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^ ! Abstract 
O 

^ I summarize and extend the qualitative results, obtained previously by inspection 

of SU(2) lattice gauge field configurations. These configurations were generated by the 
. Wilson action, then transformed to a Landau gauge and smoothed by Fourier filtering. 

This leads to sharp peaks in field strengths and related quantities, the characteris- 
J> I tics of which are very well separated from a background. These spikes are caused 



o 



by gauge singularities. Their densiis determined as 1.5/ fm , with very good scaling 



' properties as a function of the bare coupling constant. The number of spikes within a 

■ configuration vanishes when approaching the deconfinement region. Furthermore, the 

■>!:::j- ■ Landau-gauging procedure becomes unique, if the probability to find a spike is much 

I smaller than unity. The relation of the spikes to the instantons obtained by cooling 

"j^ I is investigated. Finally, a correlation between the presence of spikes and the infrared 

behaviour of the gluon propagator is demonstrated. 

Oh! 
(D . 

1 Introduction 

^ I If one transforms a single configuration of SU(2) lattice gauge fields to a Landau 

' gauge and then removes the high momentum Fourier components of the gauge 

fields by some exponential cut-off ("Fourier filtering", F.f., or "smoothing"), one 
observes the following phenomenon P]: For sufficiently large lattices and/or for 
sufficiently large bare couphng constants, at a couple of lattice positions the 
gauge fields show a rapid variation in all space-time directions and for all colour 
components. These variations lead to narrow spikes in the Wilson action density 
S{x) and in the topological charge density^ As an example, in fig. [T]the 

charge density q{x) is shown for a plane of a large lattice. The cones have been 
generated by the VRML-language [2J and visualized by the browser GLVIEWjS]. 
The positions of almost all the spikes do not vary in the different gauges which are 
due to different Gribov copies For sufficiently modest filtering, the spatial 
extension of these spikes amounts to one or two lattice units, which roughly 
corresponds to 0.05/m. 

Due to Fourier smoothing, the quantities S{x) and q{x) are gauge dependent, 
and their physical interpretation is not straightforward. The observation that, 
within the spikes, one observes self-duahty, i.e. q{x) ~ ±S{x) within a 10% 
accuracy, might suggest that we observe narrow instantons. Two facts contradict 



^The normalization of these two quantities will be given in eqs. II 31 and [HI 




Figure 1: The topological charge density on a lattice plane running through a 
typical spike. The direction of the cones indicates the orientation of the charge 
in three-dimensional colour space. The lattice size is 48^ x 64, and jS = 2.85. 
The cut-off for the Fourier transformed gauge fields is = 0.5a^. 



this interpretation. First, one observes that the distribution in terms of S{x) - 
summed over a few lattice sites around the spikes- is strongly peaked for large 
action values. Those correspond to small-sized instanton-like objects. For details 
see section El On the contrary, several lattice studies in SU(2) predict a size 
distribution with a maximum for sizes in the range from 0.2 fm [5j to 0.4 fm 
jH] and around 0.3 fm or higher in SU(3) 0. Secondly, if the configuration 
is smoothed by the standard cooling technique |H], then one indeed observes 
instantons at the position of the spikes, but their size is much larger than that 
of the spikes. 

If one transforms to a Landau gauge prior to cooling, one can directly observe 
that the gauge fields are singular (regulated by the lattice) at the positions of 
the spikes. All observations point towards the interpretation that the Landau 
gauging leads to singularities in the gauge fields, if the lattice parameters are 
within the confinement region. A gauge singularity at the origin can be created 
by imposing the gauge transformation^ 



on a relatively smooth gauge field configuration, e.g. on a large-size instanton in 
the regular gauge^. Then, close to the origin, there occurs a cancellation between 
the derivative terms in the action and the cubic and higher gauge field terms. 

This cancellation will be mutilated by filtering, which acts differently on 
the two contributions. Accordingly, a singular gauge field distribution will de- 
velop, under F.f., a sharp peak in the action and in q{x), even if both quantities 
are smooth in x without filtering'^. Thus, the plaquette action and q{x) -after 
smoothing- no longer possess a definite physical meaning. They have to be un- 
derstood as a measure for unexpected things that happen while gauging. 

In this paper, I will concentrate on the phenomenology of the spikes, i.e. on 
their density and on their correlation with other nonperturbative phenomena. 
The following topics will be treated: 

1. For large f3= 4/(7q, the spikes are very well defined objects, if the action (or 
the topological charge density) -summed over a few lattice units around the 
action maximum- is taken as a probe. This is demonstrated by the histogram 
in fig. 121 which is based on 0(30) configurations on large lattices and large /3. 
The histogram displays a well defined valley between a peak containing the 
spikes (all satisfying self-duality g(x)~ ±S{x)), and a background which 
contains distinctly different parts of the lattice configuration. This allows 
to determine the density of spikes as 



Since the average number of spikes per configuration is well determined, i.e. 
without a significant systematic error, the density of spikes in physical units 
can be used for a scaling study as a function of (3. Very good scaling with 
the string tension results ^2] is observed in the interval 2.70 < (3 < 2.85. 
For details see section El The relation of the number of spikes and the 
deconfinement transition will also be discussed in that section. 

^Thc origin x = has to be chosen apart from a lattice site, r arc the Pauh matrices. 

■^The generahzation to several singularities at different positions is obvious. 

"^This phenomenon has been checked by analytic calculations in the continuum by H. Joos. 
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2. The density of spikes is correlated with the number of Gribov copies which 
are found during the Landau gauging process. If the lattice is chosen so small 
that the probability for finding one or more spikes in the configuration is 
much less than unity, then, empirically, the gauging procedure is practically 
unique. This means that the Landau gauging algorithm ends up -with high 
probability- preferentially in the same gauge, i.e. one rarely observes the 
appearance of Gribov copies jl] if only a finite number of gauging trials is 
performed (a trial is defined by applying a random gauge prior to the Landau 
gauging). The observation is that if copies are found, the probability to find 
spikes is enhanced too. This positive correlation between the appearance of 
spikes and copies will be demonstrated in section |3] 

3. It is instructive to study alternatives to the Fourier filtering technique, no- 
tably the cooling technique |8j. In sectional I will show that there is an 
interesting connection between the outcome of cooling and of F.f. If one 
starts from a configuration in the Landau gauge, finds the position of its 
spikes and then cools the unfiltered lattice configuration, one observes the 
following: The absolute values of the gauge fields, |'U^(a;)| (see eq. El), de- 
crease under cooling almost everywhere, except at the positions of the spikes. 
There, gauge singularities show up. The local action maxima, which always 
accompagny the singularities, are broad in most cases. 

4. The fact that the spikes manifest themselves as zeros in the gauge fields -in 
all directions and for all colours- suggests that the gluon propagator and 
other correlators of the gauge fields behave differently if spikes are present 
in a configuration or if not. This is the it will be shown in section IHl 



2 Landau Gauge and Filtering 

Here a couple of topics relevant for gauging and filtering will be discussed. The 
Landau gauging algorithm on the lattice maximizes the sum of the "large" SU(2)- 
components Uq^^{x) in the link representation 

U^{x) = Uo^^{x) + iTu^{x) (3) 

with 

<^(a;) + ^'(^) = l- (4) 
Thus, one searches numerically for maxima of 

FiU)=J2uoAx)- (5) 

X,fJ, 

under local gauge transformations U. If the maximum condition is fulfilled, one 
has 

Ui,{x) -Ui,{x- Cf,) = 0, (6) 

which corresponds to 

d^A^ = (7) 

in the continuum. At the center of an instanton in the singular gauge, eq. [7| 
makes no sense. This fact requires to describe an instanton around the center by 



gauge fields in the regular gauge and at infinity by the singular gauge. On the 
lattice, however, the singular gauge is tolerable, if the center is not located at a 
lattice site, and lattice configurations in the Landau gauge seem to choose this 
option. 

The Gribov ambiguity 4j implies, that there exist -for sufficiently large lat- 
tices in physical units^- many different maxima where eq. IHlis satisfied. A priori, 
it is not clear whether gauge dependent average values of link variables, e.g. the 
gauge field propagator, depend strongly on the special Landau gauges, which are 
characterized by the different values of the gauge function eq. El for a given con- 
figuration. Empirically, the gauge field propagator does not depend strongly on 
the special gauge. This can be demonstrated by studying the correlation between 
numerical results for the gluon propagator and the values of the gauge function 
F{U) for different gauges. The latter can be generated either by different gaug- 
ing algorithms (e.g. overrelaxation or conjugate gradient) or by different random 
gauges, applied before the iterative gauging starts. The following numerical re- 
sults have been obtained: 

For a lattice of size 24^ x 32 at (3= 2.5, corresponding to a physical size 
(1.9)'^ X 2.5 /m^, 20 random gauges were used for each configuration, and the 
values for the gluon propagator at the smallest momenta were ordered with re- 
spect to the values of the gauge function, i.e. whether F{U) was smaller or larger 
than the average value (which is determined for the individual configuration). No 
significant difference between the two data sets could be observed, with an over- 
all relative accuracy of 5% for the gluon propagator at = 0, and with 1.5% 
at > 0. The upper limit for the difference between one of the ordered data 
sets and the average over the random gauges is 2% for > 0. The stochastic 
noise, due to measuring on different configurations, is comparable to that due to 
measuring at different gauges, with no significant correlation between the value 
of the gauge function and the value of the propagator. 

The Landau gauging drives the configuration to the region of small vector 
components u^{x) in the sense that link variables with uo^^{x) < become very 
rare. For presently available coupling constants, j3 < 2.85, this goal has almost 
been achieved. This allows for a linearization and subsequent Fourier transfor- 
mation of the gauge field variables. Besides the trivial linearization, 

u^{x) A^{x), (8) 

I use a stereographic transformation in order to minimize errors due to residual 
negative uo,^{x). There I define gauge fields by 

M^y'2/(l + uo,^) = w^{x) A^. (9) 

This stereographic transformation brings the "south-pole" Uq = —1 not to infin- 
ity, but to = 4. The Fourier filtering is applied to the gauge fields A^{x), (see 
eq. ITH|) and after suppression of high momentum Fourier amplitudes, the gauge 
fields are transformed back to SU(2)-variables. After an eventual smoothing of 
the fields, =^ A' ^ by F.f., normalized link variables will be reconstructed by 

^In section0]we will see, that the critical size is about (0.9/m)^. 



the obvious inversion: 



= (10) 

u'^ = Jl - wl/A w„ (11) 



u'o,^ = -sign for «5^^ > ^. (12) 

The differences between the two methods of hnearization amount to 5% for the 
field strengths, with no impact on the existence and properties of the spikes. 

With these variables, the action and the topological charge density will be 
defined in the following way: First of all, the field strength tensor F^^{x) is 
calculated by averaging the link products around the 4 plaquettes p, a with p,a ^ 
fj., u, which are connected to x. Given Fj^{x)^ electric fields E^°'\x) and magnetic 

fields B^°-\x) are defined, and the action is 

S{x) = -Y,{E^''\x)E^''\x) + B'^''\x)B^''\x)). (13) 
For convenience, I define 

q{x) = J2E^''\x)B^''\x) (14) 

a 

and call it topological charge density, inspite of the nonstandard normalization. 
For self-dual objects, one has S{x) = ±q{x). 

Due to the slight non-locality of the operator q{x), it should have a negative 
correlator [H] 

< qix)q{0) > < (15) 

only for x"^ > 4a^. This property has been successfully checked for 10 configu- 
rations on large lattices. The negativity breaks down for filtered and for cooled 
configurations^. This can be read off, for instance, from fig. 16 of ref. [H], and it 
also has been verified for the cooled configurations used in section El 



3 Density of Spikes 

Here it will be shown that the occurence of spikes is a well defined phenomenon 
in the sense that the spikes are separated from a background of action maxima 
by a deep valley, such that their density can be determined without a significant 
systematic error. This is true, especially, for large values of /?. For smaller 
values, the depth of the valley diminuishes somewhat, without a severe loss of 
significance. Thus, it can be tested quite well whether the density of spikes scales 
under a variation of /3 in accordance with the string tension, inspite of the limited 
number of spikes identified so far^. 

An interesting question is how the average number of spikes varies with the 
lattice volume, if f3 is kept fixed. The preliminary evidence is that the density 

^Obviously, the negativity excludes the dominance of four-dimensional, "sign coherent structures" in the 
true vacuum, as has been emphasized by (lOj. 

'''This is true because the change of scale enters the density with the 4*'* power. 



is roughly independent of the volume, and that the parameters of the critical 
volume, where one expects just one spike per configuration, are close to the 
deconfinement transition. This will be discussed at the end of this section. 
The following quantitative analysis is based on 

A) 34 configurations on a large lattice with size 48'^ x 64 at (3=2.85, with twisted 
boundary conditions [H], and 

B) 25 configurations on a smaller lattice with size 32^ x 48 at /5=2.70, also with 
twisted boundary conditions. 

Each of these configurations has been gauged to a Landau gauge and then 
Fourier filtered, with the following substitution for the momentum space ampli- 
tudes 

Mq) = A^iq) eM->^Y)- (16) 

For identifying spikes, I use = 0.5a^. This modest filtering already leads to a 
strong suppression of the perturbative noise of the A^. The filtered gauge fields 
were transformed back to x-space and subsequently reunitarized (see eq. [HH). 

For each of these filtered configurations, I have identified 0(500) local maxima 
of the action density S{x). Then, S{x) has been summed within a radius of 
R = 2a around the positions Xm of these maxima, 

E„(A2,i?)= J2 Six). (17) 

X 

For case A, the distribution of the values of Sm(A^,i?) is shown as the solid 
line histogram of fig. |21 The distribution is normalized by the physical volume 
V{(3) of the lattice (see below). It is remarkable how clearly the peak at large 
Em(A^, R) is separated from the background at smaller values, with a dip around 
Eto(A^, R) = 8. Furthermore, the data for case B are shown in fig. El as a broken 
line. In order to compare these data with those of case A, one has to fix the 
ratios of the physical volumes. For this, I use 

" ' a{(3 = 2.85)4 X 483 X 64 ^ ^ 

with dg a{p = 2.70)/a{p = 2.85) = 1.6. 

The fact that the dip is less pronounced at (3 = 2.7 than at the higher value 
of (5, is, most likely, not due to a statistical fiuctuation. Data at (3= 2.5 show 
that the peak for high values of Sm(A^,i?) changes in shape towards a plateau. 
At low P, the spikes are probably not as stable as at larger /?, and a snapshot of 
a single configuration will find both well developed spikes and decaying/growing 
ones. 

It is natural to define the average density of spikes per configuration by the 
number of objects with Em(A^, R) > 8. In case A, this gives 

NspikesiP = 2.85) = 6.6 ± 0.5 spikes per configuration, (19) 

and, using the string tension scale [T2j a(/3 = 2.85) = 0.028/m, 



PspikesiP = 2.85) = 1.52 ± 0.2 spikes / fm\ 



(20) 
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Figure 2: Histograms for summed Wilson action S(A^, R) for (3 = 2.85, lattice 
size = 48^ X 64 (full line), and for /3 = 2.70, lattice size = 32^ x 48 (broken 
line). The summed Wilson action E(A^,i?) and R are defined in eq. El The 
histograms are normalized to V(/?), the physical volume of the lattices. 



The corresponding results for case B are: 

NspikesiP = 2.70) = 9.6 ± 0.6 spikes per configuration (21) 

and 

Psp^kes{|3 = 2.70) = 1.51 ± 0.2 spikes//m^ (22) 

Thus, eqs. 1201 and EH show with good accuracy, that the density of spikes 
obeys scaling according to the non-asymptotic variation of the string tension, 
and that the density itself is a rather dilute one. The spikes themselves are 
better defined for large values of /3 than for small ones, and this implies, together 
with the scaling of the density, that they are not a remnant of the strong coupling 
regime, which would become irrelevant in the continuum limit. 

It remains to show how the number of spikes varies when -at fixed (3- the 
lattice size is chosen so small that one approaches the deconfinement region 
I take the variation of certain glueball masses as a signal for (de-) confinement, 
especially the mass ratio of Af and E~^, which shows a step |131 IT^ as a function 
of the lattice size L around 

= m^+L ^ 7. (23) 

^Tlic details of the transition to this region depend on the lattice geometry, and the present geometry is 
not optimally chosen for this purpose. 



Combining this with the ratio of the root of the string tension, V^, to the mass 
m^+, \/K/m^+ ^ 0.25, and with the value of the string tension ^2] a'^K(l3 = 
2.85) = 0.0040, I find for the deconfining lattice size 

L,tep = 27.7a. (24) 

Now, the important question is how many spikes have to be expected on a lattice 
of size 27.7'^ x 37. As will be shown in the next section, the density of spikes 
is (within an accuracy of 10%) independent of the volume. Thus, from eq. 
we have to expect 0.73 ± 0.15 spikes per configuration. This value allows to 
conclude that at the deconfinement transition the average number of spikes falls 
below unity. 

4 Gribov Copies and Spikes 

The existence of Gribov copies is made evident by applying random gauge trans- 
formations prior to the Landau gauging. Then, on large lattices, the final value 
of the gauge function F{U), eq. |31 will differ for each set of random numbers. It 
remains to be seen how small the lattices have to be such that the gauging pro- 
cedure is practically unique in the sense that for a large number of random trials 
one always ends up in the same gauge. Here I will show that the naively estimated 
number of spikes provides a good measure for the transition from non-unique to 
unique gauging. 

In detail, it will be studied how the appearance of Gribov copies is correlated 
with the existence of spikes. The study is restricted to a lattice size close to the 
critical one for which the expected number of spikes per configuration is just one, 
in the following sense: Given the number of spikes, Ngpikes, on a physically large 
lattice, Viarge, I define a critical volume by 

^crit ^large / -^spikes 

Thus, the number of spikes on Vcrit were one if the density of spikes were inde- 
pendent of the volume, which, of course, need not be the case a priori. 

A lattice slightly smaller than the critical one seems to be best suited to study 
the correlation between the two phenomena. The reason is that for lattices which 
are much larger than the critical one, the number of Gribov copies cannot be 
determined reliably, since for every random gauge one will end up in a different 
copy. On the other hand, for lattices much smaller than the critical one, copies 
are very rare, and it will take too many sweeps to find any. 

At P= 2.85, I found the average number of spikes on the large lattice, eq. UHl 
and accordingly we have 

V,„t = (0.90/m)^ (26) 

Thus, at /5 = 2.85 and on a lattice of size 24^ x 32, which corrsponds to a volume 
V = (0.72^)^, we would expect 0.41 spikes by naive geometrical scaling. The 
simulation, described in the following, shows that the average number of spikes 
per configuration is 

NspikesiP = 2.85, 24=^ X 32) = 0.37 ± 0.03 spikes per configuration. (27) 



(25) 



This is well compatible with the number 0.41 expected from naive geometrical 
scahng. I conclude that -within reasonable accuracy- the density of spikes is 
independent of the lattice volume. 

For the above lattice, the number of Gribov copies for a given lattice config- 
uration turns out to be strongly correlated with the average number of spikes 
observed on this configuration. The observation is based on two long sequences 
with about 200.000 updates for each sequence. The approximate number of Gri- 
bov copies has been determined after every 1.000 sweeps. For this purpose, 12 
random Landau gauges were performed, and the 12 values of the gauge function 
were compared and the numbers of different gauge functions, found in this way, 
were recorded. Simultaneously, the number of spikes was found by filtering the 
configuration with a cut-off = 0.5a^ and by searching for local action maxima, 
which exceeded the background by a factor three or more. For such maxima, 
\q{x)\ agrees with the action within an accuracy of 10 %. The results for copies 
and spikes are the following: 

• Configurations without Gribov copies follow each other in long sequencies, 
i.e. for 0(10.000) lattice updates. When spikes occur, they preferentially 
show up in all random gauges, and very often more than one spike is found. 
Thus, spikes and Gribov copies occur, during the process of updating, in 
lumps. 

• The appearance of Gribov copies is not exactly correlated with the number 
of spikes, i.e. not on a 1:1 basis for each configuration. There are configu- 
rations without copies^, but with spikes, and vice versa. 

• If both the number of copies and the number of spikes per configuration are 
smeared over a few adjacent measured configurations (just to improve the 
presentation), a striking correlation shows up. This is demonstrated in fig. 
El On the left-hand side, the first 100.000 sweeps are shown, with the slim 
curve giving the number of copies, and the fat curve giving the number of 
spikes (averaged over the 12 random gauges). The results from the next 
100.000 sweeps are plotted on the right hand side. 

It is evident, that the appearance of Gribov copies and of gauge singularities 
both do not stop abruptly when the lattice becomes small. Presumably, both 
phenomena are connected with the tunneling of configurations from one " smooth" 
state to another one. 

5 Cooling a Gauged Configuration 

The nature of the spikes can be studied further if smoothing by F.f. is confronted 
with the results of the cooling technique. By the latter, the quantum noise in 
gauge invariant observables is reduced during the approach to metastable minima 
of the action. In this way, non-perturbative aspects of the configuration, like 
instantons, may emerge. Since instanton configurations are close to local minima 
of the action, the distinction between genuine instantons and artifacts of the 
cooling procdure, is subjective. A comparison of the gauge fields etc. of a cooled 

Of course, it is possible that some copies might be found if more random gauges had been tested. 
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Figure 3: Averaged number of Gribov copies per configuration (slim points) 
and averaged number of spikes per configuration (fat points) as function of 
number of sweeps in long updating sequence. The lattice size is 24^ x 32, and 
P= 2.85. The Fourier cut-off is = 0.50^. 



configuration with the same configuration modified by F.f., may be helpful to 
understand in what sense instantons exist on the lattice. 

For this purpose it is necessary to quantify the impact of F.f. on instantons. 
This depends on the gauge. For standard instantons in the regular gauge, the 
effect of filtering is quite modest: The peak action of an artificial instanton of 
radius 7a (this corresponds to a radius p = 0.20 fm at /5= 2.85) will be reduced by 
F.f. not too strongly, if the filtering is done with a strong cut-off, i.e. = 4.0a^. 
The reduction amounts to a factor 

S{x = 0, A^ = 4:.0a^)/S{x = 0, A^ = 0.500^) ^ 0.49. (28) 

Thus, F.f. will not seriously affect those instantons, which are the most numerous 
ones in SU(2), according to ref. [S]. Smaller ones will feel a stronger reduction 
in peak height, such that they essentially look like broader instantons. 

On the other hand, as has been stated already, instantons in the singular 
gauge will show up, under F.f., as spikes and, therefore, cannot be missed. I as- 
sume, more generally, that F.f. will not reduce the number of visible instantons, 
if the creation and observation of spikes is properly taken into account. Further- 



more, F.f. will not create quasi-instantons nor improve the self-dual properties 
of maxima, in obvious contrast to the cooling technique. 

I have studied correlations between spikes and cooled objects in two ways. 
First of all, around the location of spikes in a given filtered configuration (denoted 
by JF), the cooled configuration is studied visually. Also, action maxima are 
studied, which arc not associated with spikes. Secondly, the inverse correlation 
will be investigated, i.e. maxima in the cooled configuration, called C, are selected 
and the filtered configuration T is visualized around the corresponding locations. 
In the final subsection, the probability for obtaining spikes by gauging cooled 
configurations is investigated. The number of spikes and -to a large degree- their 
positions seem to be independent on the amount of cooling. 



5.1 Spikes =^ Cooled Configuration 

The cooling sequence starts from the same Landau gauged configurations as the 
filtered ones. Four configurations have been cooled both with 30 and 100 steps, 
using "strong" coohng. There, each step rotates each link variable to the local 
action minimum. A total of 37 spikes has been found, out of an ensemble of 128 
action maxima. 

For the spikes, the resulting correlation between the filtered configuration and 
the cooled one is quite simple: At the location of the 36 spikes one finds, first of 
all, always gauge singularities in C, and, secondly, one observes that the spikes 
are, apart from 7 cases^°, associated with instantons. 

The first observation means that the cooled gauge fields show clean 

peaks around the position of the spikes, including sign fiips in all colours and 
directions, with a diameter of 3 or 4 lattice unit. Outside of the peaks, the 
gauge fields are noisy and not particularly small. This is due to the fact that 
cooling does not minimize the gauge fields, and Landau gauging still preserves 
their perturbative fiuctuations. The positions of the gauge singularity and the 
nature of the gauge field peaks are independent of the number of cooling steps, 
as has been checked by comparing the fields for 30 and 100 cooling steps. 

The second observation means that at the position of a spike, pronounced 
maxima of the action and of the topological charge density show up in C. The 
topological charge density has the opposite sign in the two cases cooling and 
filtering^^. After cooling, the hights of the action-peaks vary strongly from peak 
to peak. The explanation is that for instantons, the height of the action-maxima 
is a strongly decreasing function of the instanton size, such that large instantons 
do not show up spectacularly under visualization techniques. Obviously, selecting 
positions by spikes does not select instantons of a particular radius. 

For action-maxima observed after F.f. which are not spikes, the situation at 
the same position in C is not very clear-cut. At some maxima one may find 
an instanton-hke object, at others there is one close by, and in many cases one 
cannot observe any activity in a neighbourhood of reasonable size. 

^°An inspection of these exceptions reveals that either the cooled maximum of q{x) is quite low, or that 
cooling leads to an annihilation of a instanton-anti-instanton pair. 

^^This is in agreement with the interpretation of spikes as a mismatch of the quadratic terms and the 
higher order terms, caused by suppression of high frequency terms via filtering. 



5.2 Maxima of Cooled Configuration ^ Filtered Config- 
uration. 



The inverse correlation, i.e. between cooled action (or maxima selected 

in C on the one hand, and between structures in J-' (obtained by filtering) on 
the other hand, is more complicated to investigate than the previous one. This 
holds because there is much freedom in the selection of maxima. It is not the 
purpose of this paper to follow the elaborate filtering techniques of ref. for 
extracting the best instanton candidates among the many maxima which show up 
during the cooling process. I simply start from the 4 configurations which have 
been "strongly" cooled by 100 sweeps, select for each configuration 32 positions 
associated with the 32 highest maxima of |^(a;)|, and investigate the properties 
of JF around the positions of these maxima. The filtering is done both with 
parameters = 4a^ and = 0.5a^. 

The first observation is that one recovers ~ 80% of the spikes among the 
first 0(30) maxima. This is not trivial. The sizes of the instantons which are 
associated with spikes, vary considerably, such that a selection via the |g(x)|-peak 
height might lead to a failure in identifying these maxima within the first few 
dozens of peaks. I conclude that the spikes are tightly correlated with instantons 
which are produced by cooling and selected among the spacially less extended 
ones. Obviously, the spikes have a significance beyond the gauge dependent 
filtering technique. 

Next, it is highly interesting to investigate those locations in the filtered con- 
figuration, JF, which on the one hand are associated with peaks in C but, on 
the other hand, have no singularity in the gauge fields at this location. In a 
fraction of about 2/3 of those maxima, one encounters also a peak in J-' both in 
the action and in |g(x)|, with identical signs of q{x) in both configurations, i.e. 
(anti-)instantons in non-singular gauge in C may be associated with candidates 
for (anti-) instantons in JF. In the other 1/3 cases, there is no significant activity 
in J^. 

Now, the crucial question is whether these candidates have the correct prop- 
erties to be unambiguously identified with instantons. This has to be doubted 
for the following reasons: 

• In JF, the orientation of the colour components in colour space is not along 
the diagonal. This fact is best recognized by visual inspection, and it will 
not be specified in detail here. 

• The electric and magnetic fields in J-' are not perfectly self-dual. In order to 
be quantitative, I consider the measure of self-dual quality, Q^"'\y), defined 
by 

Q^'\y)= ^ 2 ^ ^ -.2 . a = 1,2, 3. (29) 
(?/) + 5W (y) 

where y is taken at the lattice positions around the maxima position x 
in C, with |x — ?/| < 4a. This distance amounts to half the radius of 
a typical SU(2)-instanton. A histogram of this quantity, taken between 
— 1 > Q'^°^^{y) > 1, is peaked beautifully at the limits ±1 for the cooled 
configuration. For the filtered configuration, however, Q^°'\y) shows a 



broad histogram for all a, extending down to Q^^^y) ~ 0. Thus, self-duality 
is not realized well in the case of F.f.. 

• The spatial sizes of the peaks in JF are much smaller than in C. Since the 
shape of the topological objects seems to be rather irregular, their radii are 
difficult to determine directly. A simpler way is to observe the reduction of 
the peak height as a function of the filtering parameter, when it is increased 
from = 0.5a^ to = 4a^. This reduction amounts approximately to 
a factor 0.03, in sharp contrast to the reduction by a factor 0.49 for an 
artificial lattice instanton (see eg. OHj). 

It has to be concluded that a close correspondence between cooled configura- 
tions and filtered ones exists mainly at the position of the spikes. In the cooled 
configuration, the gauge singularity -the origin of the spike- is preserved, and 
an instanton-like maximum of the action etc. has developed. Action-maxima in 
C, which are not associated with spikes, and which may therefore be displayed 
in the filtered configuration, do not reveal the typical properties of instantons. 
Of course, by a more refined search among the many maxima in C one may 
eventually find better instanton candidates. 

5.3 Gauge and Cooling (In-)Dependence of Spike Posi- 
tions 

In the previous subsection, it has been stated that almost all spikes in the filtered 
configurations are found at a position close to maxima in the cooled configura- 
tions. The relevant maxima are those with the largest values of the action density. 
This is rather important, since the positions of these cooled maxima are gauge 
independent. Thus, also the position of spikes has a gauge invariant meaning, at 
least in some probabilistic sense. 

Furthermore, in the process of cooling, most small scale fluctuations are elim- 
inated which, in principle, could induce the gauge singularities. The presence 
of such " dislocations" which can contribute to the topological charge [HI E| is 
a drawback of using the Wilson action both for updating configurations and for 
cooling. 

In the following, the effect of cooling will be studied once more in a different 
way. In determining the correlation between the cooled maxima and the spikes, 
the latter were defined by Landau gauging a configuration which had all short 
range fluctuations undamped. Here, the order of cooling and gauging will be re- 
versed, i.e. a non-gauged configuration will be cooled with up to 10 strong cooling 
steps. This eliminates all plaquette values smaller than 0.9. If one then finds 
a Landau gauge and applies filtering to this cooled configuration, one finds ap- 
proximately the same number of spikes as compared to the case without cooling. 
Most of the spikes show up at the same positions in the two cases. 

In detail, on a lattice of size 24^ x 32 at /5 = 2.65, 20 configurations have been 
cooled with 3 steps (case (a)), and 20 configurations with 10 steps (case (b)). Af- 
ter a random gauge had been applied, first the Landau gauging and then filtering 
were performed. In case (a) one finds 93 spikes, whereas for the corresponding 
uncooled configurations one finds 87 spikes. In case (b), the corresponding num- 
bers are 99 and 91. When the spike positions are compared between the cooled 



and the uncooled configurations, it turns out that in 15 configurations out of the 
20 cooled ones (in case (b)), more than 50% of the spikes show up at the same 
positions as in the uncooled case (within 1 or 2 lattice units). 

Because of the strong suppression of plaquettes with trace values < 0.9, this 
correlation of spike positions makes it rather unlikely that the spikes are induced 
by dislocations, in so far as these are characterized by large negative plaquette 
values. 

6 The Gauge Field Propagator and Spikes 

It is evident that close to a spike, the gauge fields show a rapid fiip of sign, for all 
colours, for all space time orientations and along all directions. Thus, it is natu- 
ral to expect that the gauge field correlators behave differently for configurations 
with spikes as compared to configurations without spikes. A convenient tool to 
study this effect is a measurement of the gluon propagator, evaluated separately 
for the two specimens of configurations. It is to be expected that the propagator 
in x-space falls off more rapidely for increasing spatial separation, when spikes 
are present than if none are around. In momentum space, this effect is repro- 
duced if the zero momentum propagator is reduced and the small -but non- zero- 
momentum propagator is enhanced when spikes are present, as compared to the 
spike-free configurations. 

This can be tested on lattices which are so small that -simply by geometrical 
considerations- the probability to find a spike is considerably smaller than one. 
For our standard value of (3= 2.85, a lattice of size 24^ x 32 has a density of 0.4 
spikes per configuration (see eq. I77j) . Since often there are more than one spike 
per configuration, the probability to find at least one spike is around 0.2, i.e. 
considerably smaller than 0.4. Thus for a total of 600 configurations which have 
been measured^^, one finds 480 configurations without spikes and 120 configura- 
tions with one or more spikes. The gauge field propagator has been measured 
separately for the two classes of configurations. 

The difference between the two values of the propagator, -spikes present or 
not- depends on the momentum q^. For = 0, the propagator with spikes is 
slightly smaller than without spikes: 

G{q^ = Q)no sp^kes ~ (1.19 ± 0.05)Giq^ = 0)spikes. (30) 

This difference is just significant. For > 0, the error is significantly smaller, 
and we have 

G{q^ ^ O.la-^) no spikes ~ (0.78 ± 0.03)G{q^ ^ O.la-^),^,,,,. (31) 

For > 1/a^, no significant difference between the propagators of the two 
classes could be observed. It is straightforward to transform the different be- 
haviour of the propagators to gauge field correlators in x-space. If one sums the 
gauge fields over 3 spatial dimensions and considers the correlator of this average 
along the t-direction, 

C{T) =< Y: A{t, x) Y,A{t + T,x)> / < (^ Ait, x)f >, (32) 

^^Out of 300.000 sweeps, after 500 updates Landau gauging and searching for spikes was performed. 



then the inequahties EDI and |M] imply that the correlator with spikes decreases 
faster with T that the correlator without spikes. The effect is small but signifi- 
cant. 

7 Conclusions 

Configurations of SU(2) lattice theory, when transformed to a Landau gauge, 
reveal special points with a scale invariant density of 1.5 points per /m^, where 
the gauge fields show a singularity, regulated by the lattice. This singularity 
clearly shows up either if the gauged configuration is cooled -with (9(50) cooling 
steps- or if the high momentum Fourier components are suppressed by some 
exponential cut-off. In detail, the phenomena are: 

• After cooling a gauged configuration, the gauge fields resemble those of a 
singular gauge (see eq. P), i.e. they shoot up and change sign at the special 
points. The gauge invariant action and the topological charge density signal 
the appearance of instantons of various sizes around those singularities, with 
the gauge fields in the singular gauge. Other -nonsingular- instanton-like 
objects show up visually with a density which is higer than that of the 
singular ones. This density drops quickly under prolonged cooling. The 
instantons associated with singularities are stable under cooling. 

• After removing the high momentum amplitudes by Fourier filtering, the 
singularities show up as spikes in the Wilson action and in the topologi- 
cal charge density. This is so because the removal of short range Fourier 
amplitudes destroys the cancellation between linear and quadratic terms in 
the field strengths. The gauge fields show zeroes, as a function of lattice 
positions, in all colours and directions, but the peaks are smoothed out rel- 
ative to the cooling procedure. When more and more Fourier amplitudes 
are removed, the spikes get, of course, broader. However, the positions do 
not vary essentially. 

• The positions of the spikes are not completely independent of the special 
Landau gauge, but almost so. This means that for different Gribov copies of 
a given configuration, almost all spikes appear at the same position, with a 
mismatch in the order of 10% (see .Ij). Furthermore, the spikes are strongly 
correlated with the -gauge invariant- positions of the leading action maxima 
which are generated by cooling. 

• The density of spikes is correlated with the appearance of of Gribov copies. 
This holds in the sense that on physically small lattices, where Landau 
gauging is almost unique and where the probability for the occurrence of 
spikes is smaller than unity, Gribov copies preferentially show up in the 
same configurations as spikes do. 

• The gauge field propagator is different for configurations with spikes as 
compared to the case without spikes. This shows up as a faster temporal 
decrease of zero momentum gauge field correlators. 

This holds if the scale as a function of (3 is taken from a measurement of the string tension. 



In summary, it is evident that the presence of spikes is strongly correlated with 
other nonperturbative phenomena on the lattice. In particular, the correlation 
with the gluon propagator implies that the presence of spikes is connected with 
the decrease of the propagator in x-space. The sign flip of gauge fields, which are 
associated with gauge singularities, intuitively provides a mechanism for a fast 
decorrelation of the fields. According to ref. PS], such a decorrelation can be 
responsible for deconfinement. A study of the correlation of the string tension 
with the spike density on larger lattices seems to be topic for future investigations. 



Acknowledgement 

The author is indebted to H. Joos, I. Montvay, G. Schierholz, R.L. Stuller and 
H. Wittig for useful discussions and for encouragement. The preparation of the 
SU(2)-configurations on the large lattice has been accomplished on 108 nodes of 
the T3E parallel computer at the NIC at the Forschungszentrum Jiilich. The 
author is grateful for granted computer time and support, especially to H. Attig. 
The development of the visualization tools has been carried out on a dedicated 
SNI Celsius workstation. The author is indebted to the DESY Directorium for 
providing access to this facility. 



References 

[1] F. Gutbrod, EPJdir, C8 (2001), 1, |hep-Iat/0011046| 

[2] A.L. Ames, D.R. Nadeau and J.L. Moreland, The VRML Source Book, 2nd ed., 
J. Wiley & Sons, New York 1997 

[3] H. Grahn, download of GLVIEW 4.3 from www.snafu.de/~hg 

[4] S. Gribov, Nucl. Phys. B139 (1978) 1 

[5] Th. DeGrand, A. Hasenfratz and T. Kovacs, Nucl. Phys. B520 (1998) 301 

[6] Ph. de Forcrand, M, Garcia Perez and I.O. Stamatescu, 
Nucl. Phys. B499 (1997) 409 

[7] A. Hasenfratz and C. Nieter, Phys. Lett. B439 (1998) 366 

[8] UKQCD Collaboration, D.A. Smith and M. Teper, Phys. Rev. D58 (1998) 104505, 
and references quoted therein. 

[9] E. Seller and I.O. Stamatescu, MPFPAE/PTh 10/87 

[10] I. Horvath, S.J. Dong, T. Draper, F.X.Lee, K.F. Liu, H.B.Thaker and J.B.Zhang, 
|hep-lat/0203027| Phys. Rev. D67 : 011501 (2003) 

[11] D. Daniel, A. Gonzales- Arroyo and CP. Korthals Altes, Phys. Lett. B251 (1990), 
559, and references quoted therein. 

[12] S.P. Booth, A. Hulsebos, A.C. Irving, A. McKerrel, C. Michael, P.S. Spencer and 
P.W. Stephenson, Nucl. Phys. B394 (1993), 509 

[13] C. Michael, G.A. Tickle and M.J. Teper, Phys. Lett. B207(1988), 313 

[14] P. van Baal and A.S. Kronfeld, Nucl. Phys. Proc. Suppl 9 (1989) 227 

[15] M. Gockeler, A.S. Kronfeld, M.L. Laursen, G. Schierholz and U.-J. Wiese, 
Phys.Lett. B233 (1989) 192 

[16] H.G. Dosch and Yu.A. Simonov, Phys. Lett. B 205 (1988) 339 



